Antiferromagnetic long range order in the uniform resonating valence bond state on 

square lattice 
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With extensive variational Monte Carlo simulation, we show that the uniform resonating valence 
bond state(U-RVB) on square lattice is actually antiferromagnetic long range ordered. The ordered 
moment is estimated to be m 0.17. Finite size scaling analysis on lattice up to lattice size of 50 x 50 
shows that the spin structure factor at the antiferromagnetic ordering wave vector follows perfectly 
the Sq^(.^^^) ~ + a(l/L)4 behavior, where L is linear scale of the lattice. Such a behavior is 
quite unexpected from the slave Boson mean field treatment or the Gutzwiller approximation of the 
uniform RVB state. 

PACS numbers: 75.10.Kt,71.27.+a 



The resonating valence bond(RVB) state plays 
an important role in our understanding of exotic 
physics in quantum antiferromagnet and the high Tc 
superconductors The RVB state is a quantum many 
body spin singlet constructed from coherent superposi- 
tion of different pairing patterns of the spins on the lat- 
tice. 

In practice, the most studied RVB states arc those 
generated from Gutzwiller projection of Fermionic or 
Bosonic mean field state with condensed singlet pairs. It 
is well known that the RVB states generated from these 
two routes differ in their property qualitatively. For ex- 
ample, a Bosonic RVB state can be either short ranged in 
its spin correlation when the spinon excitation spectrum 
has a full gap, or , long range ordered when the Bosonic 
spinon condense. On the other hand, a Fermionic RVB 
state can have more choices for its spin correlation. Be- 
side being short ranged with an exponentially decaying 
spin correlation, the spin can also choose to stay in vari- 
ous kind of critical states. For example, a Fermionic RVB 
state can have a large spinon Fermi surface or various 
kind of node structure in its spinon dispersion. The spin 
correlation in such critical state generally follows power 
law behavior. However, it is generally believed that it is 
difficult to describe state with long range magnetic order 
within the Fermionic RVB scheme without breaking the 
spin rotational symmetry. Such a difference between the 
Bosonic RVB state and the Fermionic RVB state stems 
from the fact that the magnetic long range order in the 
Bosonic RVB state actually originates from the conden- 
sation of the Bosonic spinon, while in the Fermionic RVB 
state, the Fermionic spinon must be paired before con- 
densation. However, a condensation of triplet pair will 
inevitably break the spin rotational symmetry. 

The relation between the Bosonic and Fermionic RVB 
state has been studied by many authors [l^, [ill- It is 
found that when the RVB ampHtude is extremely short 
ranged, or, extending only between nearest neighboring 
sites, and the lattice is planar, a Fermioinc RVB state can 
be transformed into a Bosonic RVB state with suitable 



redefinition of the RVB amplitudes. For Fermionic RVB 
state with longer range RVB amplitudes, their relation 
with Bosonic RVB state is still elusive. More generally, 
it is found that on the bipartite lattice, both Bosonic and 
Fermionic RVB state satisfy the well known Marshall sign 
rule for unfrustrated antiferromagnet, if they are both 
derived from a bipartite mean field Hamiltonian. 

An analytical treatment of the RVB state is diffi- 
cult. The most commonly adopted method to study 
the Fermionic RVB state is the slave Boson mean field 
theory. Beyond the mean field theory, there is also the 
more elaborate effective gauge field theory treatment 0- 
However, it is well known that the mean field theory 
may fail qualitatively as the no double occupancy con- 
straint is relaxed to a global one in such treatment. The 
gauge fluctuation correction, which can in principle im- 
prove the result of the mean field theory, are accounted 
for only at the Gaussian level in most cases and missed 
the singular gauge fluctuation effect at the lattice scale. 
As a variational state, the Fermionic RVB state is also 
studied by the Gutzwiller approximation [H], which try to 
estimate the effect of the local constraint by correcting 
the mean fleld expectation value with some correction 
factor. However, the Gutzwiller approximation can not 
produce any qualitatively different prediction from the 
mean field theory as to the spin correlation in the RVB 
state 0. 

The variational Monte Carlo(VMC) method is a di- 
rect way to study the properties of the Fermionic RVB 
state. It is found that in certain cases the VMC result 
can be quite different from the mean field prediction. For 
example, on bipartite lattice, a well defined sign struc- 
ture(the Marshall sign rule) may emerge out of a general 
mean field state through the Gutzwiller projection, even 
if the mean field state breaks the time reversal symme- 
try and has a complex valued wave function. Such a 
sign structure will restore the time reversal symmetry in 
the projected state and remove from the system possible 
topological degeneracy. 

In this paper, we present another example in which 
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the mean field prediction is qualitatively changed by the 
Gutzwiller projection. We show the well know uniform 
RVB state on the square lattice with a large Fermi surface 
for the spinon actually describes a state with magnetic 
long range order. 

The uniform RVB state first appears as a variational 
state for the high-Tc superconductors in Anderson's orig- 
inal paper on RVB theory of the cupratesQ- The state 
has a large nested Fermi surface for the Fermionic spinon. 
Although it is later found that the d-wave RVB state 
is energetically more favorable than the uniform RVB 
stateQ, the uniform RVB state is nevertheless still an 
interesting state of its own right. 

At the mean field level, the uniform RVB state is de- 
scribed by a filled Fermi sea with a nested Fermi surface. 
The spin correlation function decays algebraically with 
distant as The spin structure factor at the antifer- 
romagnetic wave vector is given by ^^^(Tr^Tr) — W' 
which TV is number of lattice sites I9f]. Thus at the mean 
field level, the uniform RVB state describes a state with 
no magnetic long range order. However, as the system 
presents a nested Fermi surface, the mean field state is in- 
stable with respect to infinite small perturbation toward 
antiferromagnetic ordering. Thus it is highly possible 
that the Gutzwiller projection may induce finite antifer- 
romagnetic order in the uniform RVB state, in a way that 
the spin rotational symmetry remains intact. 

The uniform RVB state is given by the Gutzwiller pro- 
jection of the ground state of the following mean field 
ansatz 

Hu-RVB (4tCj,ct + h.c), (1) 

whose mean field dispersion is given by = — 2(coskj; + 
coskj,). The uniform RVB state can then be written as 

|C/-RVB)-Pg n 4/-kj0)- (2) 

k<kF 

where kp denotes the Fermi momentum of the spinon 
and is determined by Ckp = 0. Pg denotes the Gutzwiller 
projection into the subspace of no double occupancy. The 
state can also be written in the form of condensed pairs, 

\U - RVB) = Pg ait ~ j)cl^c]^ j |0), (3) 

in which a{i — j) = X]k<kF e*'^'^*"^^. Due to the bipar- 
tite nature of the square lattice, it can shown that the 
uniform RVB state satisfy the Marshall sign rule for a 
unfrustrated antiferromagnet^, Il2l|. More specifically, 
the wave function of the uniform RVB state in a Ising ba- 
sis is real up to a global phase factor and its sign is given 
by (— , where denotes the number of down spins 
in A sublattice. 



To detect possible magnetic long range order in the 
uniform RVB state, we focus on the spin structure factor 
at the antiferromagnetic ordering wave vector 

At the mean field level, the spin structure factor can be 
found to be given by 5q=(7r,7r) — On the other hand, 
when the no double occupancy constraint is enforced by 
the Gutzwiller projection, the result is quite different. 
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FIG. 1: The spin structure factor of the uniform 
RVB state at q — (tt, tt) on a L x L square lat- 
tice as a function of the inverse linear scale of the lat- 
tice, 1/L. Presented in this figure is the results for 
L = 6,8,10,12,14,16,18,20,22,24,26,30,36,40,44 and 50. 
Periodic-antiperiodic boundary condition is adopted in the 
calculation to satisfy the closed shell condition. The error 
bars are smaller than the size of the symbols. 

In Figure 1, we show the VMC result for the spin struc- 
ture factor for lattice with size up to 50 x 50 as a function 
of the inverse linear scale of the lattice, 1/L. In our calcu- 
lation, we have used the periodic-antiperiodic boundary 
condition to satisfy the closed shell condition. We have 
used up to 4 X 10^ statistically independent samples to 
estimate the spin structure factor. Each sample is drawn 
from the Markov chain with N steps of local updates. 
The statistical error estimated from the data is smaller 
than the size of the symbols in the figure. The spin struc- 
ture factor of the uniform RVB state clearly extrapolates 
to a finite value in the thermodynamic limit. 

To see the trend of the spin structure factor more 
clearly, we fit the data with the formula 5q=(Tr,7r) = 
So + a(l/i)'^. The best fit to the data is found to be 
(3 — J and 5*0 = 0.038. Thus the ordered moment in the 
uniform RVB state is given by m = — 0.17, which is 
quite considerable as compared to the prediction of the 
spin wave theory for the Heisenberg model on square lat- 
tice. In Figure 2, a linear fit to the spin structure with 
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(l/L)* as the coordinates is shown. The hnear fit is seen 
to work extremely well from the smallest lattice we have 
studied (6 x 6) directly to the largest lattice 50 x 50. 
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FIG. 2: The spin structure factor of the uniform RVB 
state at q = (tt, tt) on a L x L square lattice as a func- 
tion of (1/L)3. Presented in this figure is the results for 
L = 6,8,10,12,14,16,18,20,22,24,26,30,36,40,44 and 50. 
Periodic-antiperiodic boundary condition is adopted in the 
calculation to satisfy the closed shell condition. The gray dot 
on the y-axis denotes the fitted value of So. The error bars 
are smaller than the size of the symbols. 

As a comparison, we also present the spin structure fac- 
tor of another well known Fcrmionic RVB state on square 
lattice, namely, the tt flux phase in Figure 3. The 7r-flux 
phase is generated by the following mean field ansatz 



H 



TT — fluX — 



= - E i^' 



'c; „Cj^a + h.c), 



(5) 



in which the phase factor is introduced to guarantee 
that each plaquette of the square lattice is threaded by a 
U{1) flux of value tt. In the 7r-flux phase, the spinon has 
a Dirac type linear dispersion and the mean field spin 
correlation also decay with distance as After the 

Gutzwillcr projection, the spin correlation of the system 
is greatly enhanced. But as is clear in Fig. 3, the 7r-flux 
phase does support an antiferromagnetic long range or- 
der. In fact, a finite size scaling analysis shows that the 
spin structure factor decays as I/L2 with the linear scale 
of the lattice and extrapolates to zero in the thermody- 
namic limit(Fig. 4). 

From these data, it is clear that uniform RVB state in- 
deed describes a state with antiferromagnetic long range 
order. This is quite unexpected from the mean field the- 
ory and it represents the first example of Fermionic RVB 
state with magnetic order. Such a order does not origi- 
nate from the condensation of spinon as in the Bosonic 
RVB state and respect the spin rotational symmetry. In 
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FIG. 3: The spin structure factor of the tt-Aux phase at q = 
(tt, tt) on a 1/ X L square lattice as a function of the inverse 
linear scale of the lattice as compared to the result of the 
uniform RVB state. Presented in this figure is the results for 
L = 6, 8, 10, 12, 14, 16, 18, 20, 22 and 24. Periodic-antiperiodic 
boundary condition is adopted in the calculation to satisfy 
the closed shell condition. The error bars are smaller than 
the size of the symbols. 
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FIG. 4: The spin structure factor of the 7r-flux phase 
at q = (tt, tt) on a L X L square lattice as a func- 
tion of 1/L^'®. Presented in this figure is the results for 
L = 6, 8, 10, 12, 14, 16, 18, 20, 22 and 24. Periodic-antiperiodic 
boundary condition is adopted in the calculation to satisfy 
the closed shell condition. The error bars are smaller than 
the size of the symbols. 



the uniform RVB state, the magnetic long range order 
can be understood as a result of the induced magnetic 
ordering in a system with nested Fermi surface. 

Following the above reasoning, it is quite interesting to 
ask the following question: does the uniform RVB state 
have antiferromagnetic long range order in other spatial 
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dimension. In one dimension, the projected Fermi sea 
state is critical and the spin correlation function decay 
with distance as indicating no magnetic long range 
order. This can be taken as the result of the strong 
quantum fluctuation in one dimensional system. As fluc- 
tuation becomes weaker in higher spatial dimension, it 
is quite possible that the uniform RVB state in three 
dimension(on cubic lattice) also have antiferromagnetic 
long range order. A study of this issue, which is numeri- 
cally quite demanding for large lattice, will be presented 
in a separated paper. 

We note a recent preprint (l3| independently confirmed 
the results of this manuscript. 
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